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ABSTRACT 

All dark matter structures appear to follow a set of universalities, such as phase-space density 
or velocity anisotropy profiles, however, the origin of these universalities remains a mystery. Any 
equilibrated dark matter structure can be fully described by two functions, namely the radial 
and the tangential velocity distribution functions (VDF), and when we will understand these two 
then we will understand all the observed universalities. Here we demonstrate that if we know 
the radial VDF, then we can derive and understand the tangential VDF. This is based on simple 
dynamical arguments about properties of collisionless systems. We use a range of controlled 
numerical simulations to demonstrate the accuracy of this result. We therefore boil the question 
of the dark matter structural properties down to understanding the radial VDF. 



Subject headings: 

1. Introduction 

A growing number of seeming universalities 
have been identified in numerical simulations of 
dark matter structures. Most of these arc in- 
tegrated quantities, such as the density profile 
(Navarro et al. 1996; Moore et al. 1998), the 
pseudo phase-space density (Taylor & Navarro 
2001), and the velocity anisotropy (Hansen & 
Moore 2006) . The cause of these universalities re- 
mains, however, essentially unknown. 

The origin of the universalities may lie in some 
fundamental property of dark matter, be it some 
statistical mechanics (Lynden-Bell 1967; Hjorth & 
Williams 2010) or optimization of some general- 
ized entropy (Plastino & Plastino 1993; Hansen et 
al. 2005; He & Kang 2010; He 2012). It may also 
be associated to dynamical effects, like radial orbit 
instability (Henrikscn 2009; Bellovary et al. 2008), 
or phase mixing or violent relaxation (Lynden-Bell 
1967; Kandrup et al. 2003). Alternatively, it could 
just be a "coincidence", since all structures have 
been built up through similar processes of merg- 
ers and accretion (Gonzalez-Casado et al. 2007; 
Salvador-Sole et al. 2007). 



A first step towards answering the question of 
the origin of the integrated universalities, is to look 
at the actual distribution of velocities. Also the 
actual shape of the velocity distribution function 
(VDF) has been suggested to be universal (Hansen 
et al. 2006), which naturally could explain all the 
integrated universalities. 

Asking the question about what dark matter 
structures fundamentally want, is different from 
asking what dark matter structures in an expand- 
ing universe actually end up doing. We will there- 
fore not be considering structures from cosmo- 
logical simulations, since their profiles often have 
merger history and environment dependent pro- 
files. We will instead consider a range of numerical 
simulations where we have better control of their 
evolution. In this way we can repeatedly perturb 
the structures in controlled manners, as well as 
giving the structures sufficiently time that phase 
mixing between individual perturbations may be 
more complete than what is the case in cosmolog- 
ical simulations. 

A non-trivial dark matter VDF also has di- 
rect implications for direct dark matter experi- 
ments (see e.g. Vergados et al. (2008); Fairbairn 
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& Schwetz (2009); Kuhlen et al. (2010) for discus- 
sions and references). 

We here present numerical evidence that the 
origin of the shape of the tangential VDF is sim- 
ple dynamics, hence supporting the idea that dark 
matter wants to follow very simple dynamical 
rules. This explains the origin of the velocity 
anisotropy profile in the inner region of dark mat- 
ter structures, with no seeming need for advanced 
statistical mechanical or generalized entropic prin- 
ciples. However, as we will point out, we are still 
left with an unknown origin of the radial VDF 
and hence also the density profile and the pseudo- 
phase space density profile are not explained yet. 

Below we will explain the surprisingly simple 
dynamical reason for the full shape of tangential 
part of the VDF, and we will perform numerical 
simulations supporting this conclusion. 

Some of these physical arguments have been 
presented previously (Hansen 2009), however, the 
simulations presented here are significantly im- 
proved. In particular we create a set of controlled 
perturbations using energy exchange reminiscent 
of violent relaxation and dynamical friction, which 
allows the particle distributions to change signifi- 
cantly, without having the structure depart from 
spherical symmetry. At the same time the struc- 
tures are analysed only after convergence to a fully 
stable configuration has been achieved. 

2. Decomposition 

Let us consider a particle moving in the smooth 
and spherical potential of many collisionless par- 
ticles. The velocity of this particular particle can 
be decomposed into three components, namely the 
radial and the two tangential components. With 
such decomposition we can consider all particles in 
a given radial bin, and get the velocity distribution 
function (VDF) in both the radial and tangential 
directions. If the structure is non-rotating, then 
the two tangential VDF's will be identical. 

It has long been known that the radial and 
the tangential VDF's are different, and physically 
this difference may seem very reasonable for equi- 
librated systems, as we will now explain. 

We will first discuss the radial VDF. Consider a 
thin spherical bin, at radius r. If we consider the 
velocity components moving outwards in the ra- 
dial direction, then those must be compensated by 



particles further out moving inwards. This com- 
pensation must depend on the particular density 
profile of the structure. This is most clearly seen in 
the Eddington inversion method (Eddington 1916; 
Binney & Tremaine 1987), from which one can eas- 
ily derive the full radial VDF from the full density 
profile. 

Now, let us instead consider the tangential ve- 
locity components. Instantaneously the compo- 
nents arc moving in the tangential plane. For 
particles with the circular speed this means that 
the component is moving in constant density and 
constant potential. Particles moving slowly in the 
tangential plane will still be near the same density 
and potential after a short time interval. That im- 
plies that the equilibrium can be achieved simply 
by having other components in the same radial 
bin moving in the opposite direction. Therefore, 
whereas the radial VDF's depend on the full radial 
density profile, then the tangential VDF's appar- 
ently do not have to concern themselves with other 
radial bins. The tangential VDF could therefore, 
in principle, be the same at all radii. 

This argument only holds instantaneously. Par- 
ticles whose tangential velocity components are 
high, will after a longer time interval, At, be 
moving outwards into lower density regions, and 
will therefore later have converted their tangential 
component into both radial and tangential veloc- 
ity components. Effectively this implies that the 
argument, that the tangential velocity component 
moves in constant potential and constant density, 
only holds for the low velocity particles. Below 
we will use numerical simulations to show that ef- 
fectively the breakdown of this argument happens 
around 70% of the escape velocity, v csc . 

2.1. Low velocity component 

We have argued above that the low velocity 
component of the tangential VDF should have 
a simple shape, which should be the shape of 
collisions-less particles moving in a constant po- 
tential and constant density. To simulate a uni- 
form medium is rather difficult using N-body sim- 
ulations, because any power or noise will induce 
gravitational collapse, which leads to a departure 
from homogeneity. Instead one can make a very 
simple analytical argument, which allows one to 
derive the VDF of a homogeneous medium. 
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Let us consider a spherical structure which has 
a known density profile, p(r). If one has /3 = 0, 
then we can use Eddington's method to derive the 
VDF, 



f(E) 



1 



/8tt 2 



d 2 p rf* 



rf* 2 y/E- * 



(1) 



where <F(r) is the relative potential as a function 
of radius, and E is the relative energy, E = * — 
mv 2 /2. Eddington's method provides the unique 
ergodic distribution function (Binney & Tremaine 
1987). We use this method to find the VDF at all 
radii. 

Imagine that the structure is particularly sim- 
ple, namely with a constant density slope, 7 = 
dlogp / dlogr , over a very large radial range. To be 
concrete, say that 7 = — 2 over 20 orders of magni- 
tude in radius, and then truncated abruptly inside 
and outside this range. Eddington's method shows 
us that this (isothermal sphere) has a VDF which 
is a Gaussian at all radii (except for the details 
arising from the truncation). 

Now, consider a more shallow slope, e.g. 7 = 
— 1.2 or —0.4. For any value of 7 we can use Ed- 
dington's method to calculate the VDF, and at 
any radius it will have exactly the same shape as 
function of v r /a r . For each density slope we use 
this method to find the unique distribution func- 
tion. 

Finally we can extrapolate this approach to 7 = 
0, which is identical to the case of constant density 
and constant potential. The shape of the VDF for 
this case is 

-5/2 
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7 = is the condition for a particle moving in 
constant density and potential, and this is there- 
fore the shape for the tangential VDF. In the case 
with 7 = there is no difference between the ra- 
dial and tangential directions, and the assumption 
of j3 — in the derivation is therefore correct, and 
the technique is thus self consistent. 

We now have an analytical expression for the 
shape of the tangential VDF at small velocities, 
and only the normalizations are unknown. These 
are the overall normalization (which must be re- 
lated to the density at that radius), and the nor- 
malization of the velocity (which is related to the 
tangential velocity dispersion). 



2.2. High velocity component 

The high velocity components are possibly even 
simpler to describe. If a velocity component is 
purely tangential at a given time, then shortly 
later it will be a combination of tangential and 
radial (unless it happens to have exactly the cir- 
cular speed). We should therefore expect that the 
shape of the tangential and radial components are 
similar at high velocities. 

There is only one complication, namely the nor- 
malization. The overall normalization must be 
identical between the radial and tangential compo- 
nents (since this is just the local density), however, 
as opposed to the low velocity component dis- 
cussed above, the normalization of the high veloci- 
ties components must be absolute, i.e. w ra d ~ v tan - 

2.3. The transition 

The transition from the low to the high velocity 
components in a real system must be smooth, and 
is probably rather non-trivial, however, for sim- 
plicity we will here make the approximation that 
the transition is abrupt, and we make no attempt 
to make it smooth. Practically, we simply assume 
that the transition always happens near 0.7w osc . 
We will address this issue further in the discus- 
sion section. 

3. Simulations 

The first simulation is a cold collapse, where 
the inclusion of substructure breaks the spherical 
symmetry. We distributed 5x 10 5 particles accord- 
ing to a Hernquist density profile (Hernquist 1990) 
with scale radius 1, and a cutoff at 200. In addi- 
tion 5 x 10 5 particles, with the same mass as the 
main halo particles, were distributed in 24 identi- 
cal subhaloes, also having Hernquist density pro- 
files, but with a scale radius of 0.5 and a cutoff 
radius of 5. The centers of the subhaloes were 
sampled in the same way as the particles in the 
main halo. The velocities of all the particles were 
initially zero. The total mass in the simulation 
was 1. We ran the simulation for 200 time units, 
which corresponds to 200 dynamical times at the 
scale radius for the initial structure. Such a cold 
collapse is similar to the simulations by van Al- 
bada (1982). 

For all the non-cosmological simulations dis- 
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cussed here, we used the parallel N-body simu- 
lation code, Gadget2 (Springel 2005). For further 
details on the cold collapse see Sparre & Hansen 
(2012). 



8000 






6000 


^^^^^^^^^^^^^^^ 




4000 






2000 














-2-10 2 
velocity 



Fig. 1. — Radial (blue stars) and tangential (green 
diamonds) VDF for three radial bins in the cold 
collapse simulation. From top to bottom the ra- 
dial bins are at 7 = —1.6, —2.0 and —2.4, and the 
bins are shifted vertically to improve readability. 
The red (solid) lines are all of the same theoreti- 
cal shape in eq. (2), which are seen to provide an 
acceptable fit in the low velocity region. 

To sample the VDF's we distribute the particles 
in radial bins with the same number of particles 
in each bin. In Figure 1 we show both the radial 
(blue stars) and tangential (green diamonds) VDF 
for three radial bins, chosen near a slope of 7 = 
-1.6,-2.0 and -2.4 (from top to bottom). We 
also show the predicted shape of the tangential 
VDF (solid line), which is clearly seen to provide 
an acceptable fit in the low velocity region. It is 
also clear that the radial and the tangential VDF's 
are very different in the low velocity region. 



To see the details of how the radial and tangen- 
tial VDF's start agreeing in the high velocity re- 
gion we plot the VDF's from the same three radial 
bins in a lin-log space in Figure 2. It is clear that 
for high velocity particles the radial VDF's (blue 
stars) are in good agreement with the tangential 
VDF's (green diamonds). Interestingly, the tan- 
gential VDF's can be approximated with the theo- 
retical solid curve for low energies, and with the ra- 
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Fig. 2. — Same as figure 1, only linear-log to see 
details of the high-energy tail. Radial (blue stars) 
and tangential (green diamonds) VDF for three ra- 
dial bins in the cold collapse simulation. From top 
to bottom the radial bins are at 7 = —1.6,-2.0 
and —2.4, and the bins are shifted vertically to 
improve readability. The red (solid) lines are of 
the theoretical shape for the low velocity region. 
For the high velocities it is clearly seen that the ra- 
dial and the tangential VDF's approach each other 
rapidly. 

dial VDF's for high velocities. And for the radially 
innermost bins (at slopes shallower than —2) this 
transition happens to be just around v = 0.7v csc . 
Only for radial bins further out than 7 = —2, it 
seems needed to use a more smooth transition. 

3.1. G-perturbations 

In order to test further the theoretical claims for 
the tangential VDF, we wish to construct a per- 
turbation/equilibration scheme, which allows the 
VDF's to change significantly, without having the 
structure depart from spherical symmetry. 

We set up structures in perfect equilibrium. 
These structures may have any density profile, 
and have zero anisotropy or follow and Osipkov- 
Merritt beta profile (Binney & Tremaine 1987). 
Now we increase the value of the gravitational con- 
stant by 20%. This increases the potential and 
makes the structure contract, and after a few dy- 
namical times a new equilibrium is reached. Next 
we repeatedly increase or decrease the gravita- 
tional constant, and between each change we allow 
the structure to phase-mix and find a new equilib- 



4 



rium. After 20 such perturbations we use the stan- 
dard value of G, and let the structure relax com- 
pletely. For further details, see Sparre & Hansen 
(2012). 

In Figure 3 we show three radial bins from 
a simulation which initially was constructed as 
an isotropic Hernquist profile. This in particular 
means that the initial conditions (before the G- 
perturbations were executed) had identical radial 
and tangential VDF's. Now, after the perturba- 
tions and subsequent relaxation there is a large 
difference between the radial and tangential VDF 
for small velocities. Instead, at high velocities the 
tangential and radial VDF quickly approach one 
another. As is also visible from the figure, the 
tangential VDF is well fitted by the theoretical 
prediction for small velocities. 
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Fig. 3. — Radial and tangential VDF for three ra- 
dial bins in the G-perturbation simulation. From 
top to bottom the radial bins are at 7 = —1.8, —2.3 
and —3, and the bins are shifted vertically to im- 
prove readability. All the red (solid) lines are of 
the same theoretical shape in eq. (2) for the low 
velocity region. For the high velocities it is clearly 
seen that the radial and the tangential VDF's are 
very similar. 

3.2. Explicit energy exchange 

Collisionless particles experience different kinds 
of energy exchange between each other, in partic- 
ular through violent relaxation (where the chang- 
ing potential implies that the particle energies 
change), and through dynamical friction (which 
transfers energy from the fast to the slower parti- 



cles). We therefore consider a perturbation where 
the spherical symmetry is again conserved, how- 
ever, we allow the particles to exchange energy 
amongst each other. This is done in such a way, 
that each radial bin conserves energy, whereby 
both density and dispersion profiles are unaffected 
by the perturbation itself. This energy exchange 
is instantaneous and the subsequent evolution is 
with normal collisionless dynamics. After each 
perturbation we again allow for sufficient phase 
mixing (see Hansen, Juncher & Sparre 2010, for 
details). After sufficiently many perturbations 
(typically 20 or 30) the structures have converged 
to a stable state, which will not change when ex- 
posed to further similar perturbations (Hansen, 
Juncher & Sparre 2010; Barber et al. 2012). 

In Figure 4 we present the VDF from three ra- 
dial bins in the final structure, which are taken at 
density slopes of 7 = —1.7, —2.4 and —3.0. Again 
we see that the final VDF's agree well between 
the radial and tangential for high velocity compo- 
nents, and that the low velocity components are 
well fitted by the analytical expression. 
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Fig. 4. — Radial and tangential VDF for three ra- 
dial bins after perturbations by explicit energy ex- 
change. From top to bottom the radial bins are at 
7 = —1.7, —2.4 and —3.0, and the bins are shifted 
vertically to improve readability. The red (solid) 
lines are of the theoretical shape for the low veloc- 
ity region. For the high velocities it is clearly seen 
that the radial and the tangential VDF's rapidly 
approach each other. 
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4. Discussions 



radii. 



We have demonstrated that three extremely dif- 
ferent artificial and controlled perturbations all 
lead to a tangential VDF, which is in good agree- 
ment with the theoretical prediction. This result 
is in good agreement with earlier studies including 
head-on collisions and galaxy formation (Hansen 
et al. 2006; Hansen 2009), and provides very strong 
evidence that the origin of the shape of the dark 
matter tangential VDF is indeed as simple as ex- 
plained in section 2.1. An important difference 
from earlier studies is that we have here been in- 
vestigating structures which have been exposed 
to controlled perturbations, and analysed only af- 
ter convergence to a stable configuration has been 
achieved. 

Let us remind the idea behind this paper. If 
we know the radial and the tangential VDF, then 
we know everything else, such as the phase-space 
density, the velocity anisotropy and the density 
profiles. For some idealized structures it is possible 
to derive the radial VDF directly from the density 
profile, and our results here imply that in that 
case we can derive the tangential VDF. That is, if 
we have the radial VDF, then we can derive the 
tangential VDF. 

In the derivation of the tangential VDF we 
made no assumptions about the anisotropy of the 
system, and the tangential VDF should therefore 
be the same for all systems and at all radii, irre- 
spective of their anisotropy profiles. This state- 
ment naturally only holds for realistic systems 
which have been perturbed and allowed to relax. 
One can always create systems in a quasi equilib- 
rium states which may even have highly different 
distribution function. Systems created away from 
an equilibrium state most often also have very 
different tangential distribution functions. How- 
ever, as we are demonstrating in this paper, all 
systems which are exposed to sufficient perturba- 
tions and subsequently allowed to relax to a quasi 
equilibrium state, will indeed have a tangential 
VDF of exactly this shape. In figure 5 we present 
the anisotropy profiles for the 3 systems consid- 
ered. In particular one sees that the anisotropy 
has a strong radial variation, going from essen- 
tially isotropic in the central region, to radially 
dominated orbits in the outer regions. And yet 
the shape of the tangential VDF is the same at all 
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Fig. 5. — The radial variation of j3 as function 
of radius for the structures considered in figures 
1-4 (solid line HJS perturbation, dashed line G- 
perturbation, dot-dashed line infall simulation). 
The structures all have a strong variation, going 
from essentially isotropic towards the central re- 
gion, to radial orbits in the outer regions. 

It has previously been suggested that pos- 
sibly the radial VDF is sufficiently close to 
a rescaled VDF resulting from the Eddington 
method (Hansen 2009). We have tested this sug- 
gestion by fitting the density profile and then using 
the Eddington method to extract the radial VDF 
at all radii. However, the resulting VDF is not an 
accurate representation of the actual radial VDF. 
This means that we are still not at the point of 
understanding the radial VDF. 

One interesting aspect of the arguments pre- 
sented for the shape of the tangential VDF is, 
that for 7 = it should be identical to the ra- 
dial one. That trend is already clear from the 
figures, namely that the tangential VDF is sug- 
gestively close to the radial one for the innermost 
bins (upper curves in all figures). To test this fur- 
ther, we selected a radial bin in the inner region 
(outside of 5 times the softening length for each 
structure, and also outside a further 30.000 par- 
ticles). The result is seen in Figure 6, where the 
tangential (dashed) and radial (solid) VDF's are 
seen to be very similar. The local density slope 
at these bins were ranging from 7 = —0.4 to — 1 
(bottom to top on Figure 6). 

We have discussed that the transition between 
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Fig. 6. — Radial (solid) and tangential (dashed) 
VDF for an inner bin (but not including approx- 
imately 35.000 of the most central particles). It 
is clear, that for all the structures the radial and 
tangential VDF are very similar. We add the 
same theoretical shape of eq. (2) in red (dot- 
dashed) lines, just to demonstrate the agreement. 
The structures included here are the cold collapse, 
two different G-perturbations and three different 
HJS2010 perturbations, covering a range of initial 
density and anisotropy profiles, in particular do we 
here present results for structures which were ini- 
tially set up with a shallow central density profile. 
Velocities are scaled by the escape velocity and 
shifted vertically to improve readability. For this 
figure we selected structures which were created 
with zero inner slope, ^initial = 0, before pertur- 
bations were applied. 

high and low velocity may be approximated as 
a rather sharp transition. This is certainly a 
good approximation for the inner region (inside 
a slope of 7 = —2). However, at larger radii it is 
clear, that a more smooth transition would pro- 
vide a more accurate representation of the tan- 
gential VDF. For all the radial bins and for all the 
structures considered in this paper, we have esti- 
mated the best velocity for the transition, t>t ra ns 5 
and it appears that this transition is in the range 
0.6vosc < ftrans < 0.8u osc . Thus, in order to avoid 
unnecessary fitting parameters, it is a rather good 
approximation to fix this at Ktrans = 0.7w csc . 



5. Conclusions 

We have demonstrated that half 'of the distribu- 
tion function (specifically, the tangential velocity 
distribution function) for dark matter structures 
can be understood from simple dynamical argu- 
ments. This implies that when we will eventually 
be able to derive the other half (namely the ra- 
dial velocity distribution function), then we under- 
stand all the properties of dark matter structures, 
including the seeming universalities of the density, 
phase-space density and velocity anisotropy pro- 
files. 

We saw that the derivation of the tangential 
VDF did not require any reference to statistical 
mechanics or generalized entropy, but instead ap- 
pears as a result of very simple dynamics. It now 
remains to derive the radial VDF, and it will be 
interesting to see if this will also be possible based 
on similar basic dynamical arguments. 
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